CORNEAL TOPOGRAPHY ANALYSIS SYSTEM 
BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a corneal 
5 topography analysis system and a method of 
performing corneal topography analysis. 

2. Description of the Related Art 

Corneal topography analysis systems are 
known which analyze the three-dimensional 

10 topography of the cornea to be examined by 
projecting a number of Placido rings onto the 
cornea, taking the Placido ring image created 
by the convex surface of the cornea by means 
of an image sensor, and finding data for the 

15 curvature of the cornea based on the taken 
Placido ring image. The corneal topography is 
displayed in terms of various kinds of color 
maps. The information about the corneal 
topography obtained by such an analysis system 

20 is utilized for planning an operation for 
correcting the cornea, for planning cataract 
surgery, and for postoperative evaluations. In 
addition, the topography is used for early 
discovery and diagnosis of keratoconus that is 

25 one of several corneal diseases. 
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As a diagnosis of keratoconus, a method of 
judging keratoconus topography using a neural 
network approach is disclosed by Michael K. 
Smolek et al. in "Current Keratoconus Detection 
5 Methods Compared With a Neural Network Approach, 
Investigative Ophthalmology & Visual Science, 
October 1997, Vol. 38, No. 11, pp. 2290-2299" 
(Reference 1), which is incorporated herein by 
reference. In this method, the presence of 

10 keratoconus (KC) and keratoconus suspects (KCS) 
are judged using a neural network and 10 indexes 
characterizing the corneal topography, i.e., 
Differential Sector Index (DSI) , Opposite Sector 
Index (OSI) , C e n t e r - S u r r o u n d Index (CSI) , 

15 Analyzed Area (AA), Cylinder (CYL), Irregular 
Astigmatism Index (IAI) , the steep axis of 
simulated keratometry (SKI), Surface Regularity 
Index (SRI), Surface Asymmetry Index(SAI) and 
the Standard Deviation of corneal Power (SDP) . 

20 Conventionally, however, only keratoconus 
topographies have been judged. The categories 
of corneal topographies and the display of the 
results of the analysis leave room for further 
i mpr o veme n t . 

25 The spacing and number of Placido rings are 
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different among manufacturers of corneal 
topography systems for calculating the indexes 
that are necessary for classification of 
keratoconus topographies. Furthermore, the 
5 number of data items obtained by edge detection 
of ring images and the data structure are 
different among the manufacturers of the systems. 
Generally, therefore, there is no data 
comparability among the manufacturers. It 

10 follows that corneal topography is analyzed 
based on point data for the data structure for 
each individual manufacturer. In addition, even 
for the same manufacturer, if the Placido ring 
structure is varied, it is necessary to change 

15 the analysis program and display program. 

SUMMARY OF THE INVENTION 
It is an object of the present invention 
to provide a corneal topography analysis system, 
20 and a method of performing corneal topography 
analysis, that improves on the classifications 
of corneal topographies, such as keratoconus 
topographies, and that can diagnose them more 
precisely. It is another object to provide a 
25 corneal topography analysis system, and method 
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of performing corneal topography analysis, that 
absorbs differences in Placido ring structure 
or data structure that would normally occur among 
different corneal topography analysis systems 
5 and which needs no amendment in subsequent 
analysis processing, display, and other 
processing . 

To achieve these objectives, the present 
invention provides a corneal topography analysis 

10 system comprising an input unit, preferably one 
that can receive input from different corneal 
curvature measuring devices, for entering 
corneal curvature data, such as corneal curvature 
data obtained by projecting Placido rings onto 

15 a cornea and taking a Placido ring image created 
by the convex surface of the cornea, or by other 
methods for measuring corneal curvature. The 
system further comprises an analysis unit for 
determining plural indexes characterizing 

20 topography of the cornea based on the input 
corneal curvature data, for judging corneal 
topography from features inherent in 

predetermined classifications of corneal 
topography using the determined indexes and a 

25 neural network, for judging at least one of normal 
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cornea, myopic refractive surgery, hyperopic 
refractive surgery, corneal astigmatism, 

penetrating keratoplasty, keratoconus, 

keratoconus suspect, pellucid marginal 

5 degeneration, or other classifications, and for 
judging its probability and grading the severity 
of keratoconus using the determined indexes. 
In addition, the system may comprise a display 
unit for displaying results of the judgments made 

10 by the analysis unit. 

As an example of one particular application, 
the analysis unit judges keratoconus cases from 
similarity to keratoconus and from severity of 
keratoconus using the determined indexes and the 

15 neural network. 

The invention also provides a corneal 
topography analysis system comprising an input 
unit, preferably one that can receive input from 
different corneal curvature measuring devices, 

20 for entering corneal curvature data, such as data 
obtained by projecting Placido rings onto a 
cornea and detecting a Placido ring image created 
by a convex surface of the cornea, or by other 
methods for measuring corneal curvature. This 

25 system further comprises a computational unit 
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for converting the entered corneal curvature data 
into a denser first data matrix by interpolation, 
removing high-frequency components from 

resulting data by frequency analysis, and 
5 converting produced data into corneal curvature 
data in the form of a given second data matrix. 

Additionally, this system comprises an analysis 
unit for judging classifications of corneal 
topographies based on the converted corneal 

10 curvature data and a display unit for displaying 
results of the judgments. 

Other objectives and features of the 
invention, including methods of carrying out the 
invention, will become evident in the course of 

15 the description thereof, which follows. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The advantages, nature and various 
additional features of the invention will appear 
20 more fully upon consideration of the illustrative 
embodiment of the invention which is 
schematically set forth in the accompanying 
drawings . 

Fig. 1 is a schematic diagram of a corneal 
25 topography analysis system according to the 
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present invention. 

Fig. 2 is a flowchart illustrating analysis 
of corneal topography. 

Fig. 3 shows graphs of data for corneal 
5 curvatures before and after noise removal and 
smoothing processing. 

Fig. 4 is a diagram illustrating a method 
of calculating the keratoconus index (KCI) . 

Fig. 5 is a block diagram of a neural network 
10 for judging the classifications of corneal 
topography . 

Fig. 6 represents an example of display of 
the results of analysis of classifications of 
corneal topographies. 
15 Fig. 7 represents an example of a display 

of a color map showing the distribution of corneal 
curvatures . 

Fig. 8 represents an example of polynomial 
expressions of a neural network. 

2 0 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Preferred embodiments of the present 
invention are hereinafter described with 
reference to the drawings. Figure 1 is a schematic 
25 diagram of a corneal topography analysis system 
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according to the present invention. This system 
includes a measuring optical system 10 for 
measuring the corneal topography. The optical 
system 10 includes a Placido disk 12 on which 
5 a large number of concentric Placido rings are 
formed about an optical axis LI, an illumination 
light source 13 for substantially uniformly 
illuminating the ring pattern on the Placido disk 
12, an objective lens 14, and a CCD camera 15. 

10 The objective lens 14 and CCD camera 15 are used 
to take the ring pattern image projected onto 
a cornea to be examined. The structure of the 
Placido disk 12 may be of the cone type. The 
objective lens 14 and CCD camera 15 act also as 

15 a viewing optical system for observing the 
anterior eye part of the eye examined. The 
topography system is fitted with a fixation 
optical system and an alignment optical system 
(not shown) . Detailed descriptions of the 

20 fixation optical system and the alignment optical 
system are omitted, since such descriptions are 
not necessary for one skilled in the art to fully 
understand the invention. 

While the measuring optical system 10 

25 described above employs a Placido disk for 
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measuring the corneal topography, the invention 
is not limited in this respect, and other corneal 
measuring devices may be used. 

The image taken by the CCD camera 15 is 
5 accepted into a video capture device 22, which 
is connected with an a r i t h m e t i c - a n d - c o n t r o 1 unit 
20 via a bus 23. The a r i t h me t i c - a n d - c o n t r o 1 unit 
20 analyzes or otherwise processes the Placido 
ring image taken by the CCD camera 15. An image 

10 control unit 21 is connected with a liquid crystal 
display 24 that can provide a color display to 
control the image displayed on the display 24 
or that can display the results of the analysis. 
A memory 25 for storing the Placido ring image 

15 and the results of the analysis, a hard disk (HDD) 
26 loaded with an analysis program or the like, 
a serial I/O 28 connected to a keyboard 29 and 
a mouse 30, a parallel I/O 31 connected to a 
printer 32, a control switch 34 having various 

20 switches, a communication port 33, and so on are 
connected to the bus 23. An external computer 
40 can be connected to the communication port 
33 through which data can be sent and received. 

A method of taking the Placido ring image 

25 and subsequent analytical processing are 
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described n ex t . During measurement, the 
measuring optical system 10 is appropriately 
aligned with the eye to be examined- Subsequently, 
a trigger switch is depressed to light the light 
5 source 13. A Placido ring pattern created by the 
Placido disk 12 is projected onto the cornea of 
the eye being examined. The Placido ring pattern 
projected onto the examined eye is taken by the 
CCD camera 15. The resulting Placido ring image 

10 is stored in the memory 25. The 

a r i t hme t i c - a n d - c o n t r o 1 unit 20 detects the edges 
of the boundaries between black and white circles 
of the ring pattern from the image stored in the 
memory 25. The results are sent to the image 

15 control unit 21. The image control unit 21 creates 
and displays a Placido ring image at the 
corresponding edge positions of the taken ring 
image, the displayed image being rimmed with an 
easily recognizable color, which is displayed 

20 on the display 24. 

If the image derived by detection of the 
edge positions and displayed on the display 24 
is acceptable for accurate measurement, the 
examiner clicks on an analysis start key 

25 displayed on the viewing screen of the display 
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24 using the mouse 30. Thus, an analysis is 
started. An analysis start signal causes the 
a r i t hme t i c - a n d - c o n t r o 1 unit 20 to execute an 
analysis of corneal topography such as the one 
5 represented by the flowchart of Fig. 2 and 
described next. 

The a r i t hme t i c - a n d - c o n t r o 1 unit 20 detects 
the edges of the boundaries between black and 
white circles of each Placido ring image as 

10 described above (step 1) . The 

arithmetic-and-control unit 20 calculates the 
corneal curvature (axial curvature) at each given 
angle based on the distance of each edge from 
the corneal center (step 2) . Details of the method 

15 of calculating the corneal curvature are 
described in J P - A - H 7 - 1 2 4 1 1 3 and U.S. Patent No. 
5,500,697. For example, assume that there are 
23 rings to be analyzed and that sampling is done 
at every one degree of the measuring meridian 

20 angle. The corneal curvature of a 2 3 x 3 6 0 polar 
coordinate data matrix is obtained. The obtained 
data for the corneal topography depends on the 
structures intrinsic to the system and on the 
method of processing, such as the number of 

25 Placido rings, the ring spacing, and the sampling 



- 11 - 



angle. The axial curvature that is one kind of 
fundamental data indicating the corneal 
topography contains a slight amount of detection 
error (or noise) in the stage of detection of 
5 the edges of black and white circles of the 
Placido ring image. 

To eliminate the dependence on the 
structures intrinsic to the system, the following 
steps are performed. First, the data for the 

10 corneal curvature in the form of the polar 
coordinate data matrix entered (or obtained) as 
described above is converted into denser 
orthogonal coordinate data ( x y - c o o r d i n a t e data) 
(step 3) . At this time, the value of the corneal 

15 curvature at each coordinate is calculated from 
the corneal curvatures at polar coordinates by 
interpolation. Various methods can be adopted 
for the in t e rpola t i ve calculations of corneal 
curvatures. For example, the space between two 

20 edge positions can be simply linearly 
interpolated. For example, if 23 black and white 
edge rings exist in an area having a diameter 
of 12 mm to be subjected to corneal topography 
analysis (edge spacing equals about 0.52 mm) and 

25 form a data structure intrinsic to the system, 
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the data is calculated by interpolation within 
the same area on the cornea and converted into 
256 x 256 orthogonal coordinates. Assuming that 
the diameter of the area on the cornea is 12 mm, 
5 a 2 5 6 x 2 5 6 data matrix is produced at a spacing 
of 0.047 mm on the cornea. 

Then, the orthogonal coordinate data matrix 
is subjected to two-dimensional fast Fourier 
transform (2D F FT ) and converted into frequency 

10 domain data (step 4) . After removing certain 
high-frequency components (step 5), an inverse 
Fourier transform (IFT) is performed (step 6) . 
This results in smoothed corneal curvature data 
from which high-frequency component noises have 

15 been removed. Fig. 3a is a graph of the corneal 
curvature data for the horizontal (X-axis) 
coordinates in the stage of interpolation 
calculation after conversion into the orthogonal 
coordinate data in step 3. In contrast, Fig. 3b 

20 is a graph of corneal curvature data obtained 
after the execution of the inverse Fourier 
transform in step 6. It can be seen that noise 
removal and smoothing have been performed. The 
frequency analysis is not limited to F FT . Other 

25 frequency analysis techniques such as wavelet 
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transform techniques may be used. 

The inverse F o u r i e r - t r a n s f o r m e d orthogonal 
coordinate data are reconverted into a polar 
coordinate data matrix (step 7) . This polar 
5 coordinate data matrix may be converted into a 
given data structure convenient for calculations 
of indexes characterizing the corneal topography 
and for map representation. For example, the data 
structure is returned to the coordinates of the 

10 edge positions of 23 rings and measuring 
meridians spaced from each other by 1° such that 
the existing analysis and display routines about 
map representations can be applied. 

Alternatively, a new data structure in which 

15 measuring meridians spaced from each other by 
1° and a given spacing of about 0.1 to 0.2 mm 
are used may be created. In this case, data arrays 
are regularly spaced and so data calculation 
about the area on the cornea is facilitated. For 

20 instance, where corneal powers on strong and weak 
principal meridians at radial positions of 3 mm, 
5 mm, and 7 mm on the cornea are found and displayed, 
it is only necessary that the corresponding data 
be specified and a calculation be performed. 

25 Furthermore, the memory capacity can be saved 



- 14 - 



compared with the case where data decompressed 
into 2 5 6 x 2 5 6 orthogonal coordinate data items 
are directly stored. Moreover, the time required 
for subsequent data computation can be shortened. 
5 By adding the steps 3 to 7 described thus 

far to the related art method, a data structure 
not dependent on the structures and processing 
intrinsic to the system such as the number of 
Placido rings can be achieved. Calculations of 

10 various indexes regarding the corneal topography 
as described below and calculations for display 
of the results of the analysis are made possible 
by a common routine. 

Analysis regarding judgment of the 

15 classifications of corneal topographies is 
described next. On the basis of the corneal 
curvature data obtained by step 7 described above, 
several indexes are calculated. For example, 
the following 18 indexes are calculated: SIMK1 

20 (Simulated Kl), SIMK2 (Simulated K2), MINK 
(Minimum Keratometry Value), SRI (Surface 
Regularity Index) , SRC (Area compensated Surface 
Regularity Index), OSI (Opposite Sector Index), 
DSI (Differential Sector Index), CSI 

25 ( C e n t e r / S u r r o u n d Index), KPI (Keratoconus 
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Prediction Index), CYL (Simulated Ke r a t ome t r i c 
Cylinder), IAI (Irregular Astigmatism Index), 
ACP (Average Corneal Power) , FAA(Analyzed Area) , 
CEI (Corneal Eccentricity Index), KCI 

5 (Keratoconus Index), CVP (Coefficient of 
Variation of Corneal Power), SDP (Standard 
Deviation of Corneal Power), SAI (Surface 
Asymmetry Index) (step 8) . 

SIMK1 and SIMK2 are powers on strong and 

10 weak principal meridians at a radial position 
of about 3 mm on the cornea. Numerical values 
higher than the normal value predict keratoconus, 
c o r n e a - t r a n s p 1 a n t e d eyes, and natively steep 
corneas. Numerical values lower than the normal 

15 value predict m y o p i a - c o r r e c t e d corneas and 
natively flat corneas. 

MINK is the lowest power of the powers at 
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cataract surgery. 

SRI indicates a local variation in the center 
of the cornea and is a value associated with PVA 
(potential visual acuity) . An increase in SRI 
means that the corneal surface within the pupil 
radius is irregular. Dry eyes, corneal 
deformation due to wearing of contact lenses, 
impaired eyes, and c o r n e a - t r a n s p 1 a n t e d eyes 
indicate high values. SRC is an index obtained 
by correcting SRI by the area of the subject to 
be analyzed. SRI is calculated by the following 
formula . 



where i the ring number, j is the position on 
a meridian, A is a constant, B is a constant, 
C is a constant, and N is the number of measuring 
points . 

ACP is the average corneal power within the 
entrance pupil. This can be regarded as the SE 
value of the cornea. If the keratometry value 
is abnormal, ACP also takes an abnormal value. 

CEI means eccentricity indicating asymmetry 



SRI = 




-C 



N 
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of the cornea, A positive ( + ) value indicates 
a prolate cornea. 0 (a nil value) indicates a 
spherical cornea. A (-) negative value indicates 
an oblate cornea. A keratoconus shows a positive 
5 value much higher than those of normal corneas. 
Corneal deformation due to wearing of contact 
lenses (CL) and myopia corrected eyes show 
negative (-) values. CEI is calculated by the 
following formula. 

10 CEI = ^~ 

N 

where i is the ring number, e± is the eccentricity 
calculated from the average distance from the 
center of the cornea at the ith ring and from 
the h e i g h t ( d i s t a n c e from the corneal vertex plane 
15 to the ith ring), and N is the number of measuring 
points . 

CVP is a numerical value indicative of an 
index created from the power distribution over 
the whole region on the cornea and calculated 
20 from the following formula. 

CVP = 1 000 x (standard deviation SD of powers 
over the whole region on c o r n e a / a v e r a g e power 
over the whole region on cornea) 

Where the value of the CVP is high, the cornea 
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shows a multifocal nature. For example, this is 
seen on keratoconi of moderate to heavy degree. 
Furthermore, it is seen during or after a cornea 
transplant surgery. 
5 With respect to OSI, DSI, CSI, IAI, FAA, 

SAI, and KPI, they are described by Maeda N., 
Klyce SD, Smolek MK, and Thomson HW in "Automated 
Keratoconus Screening With Corneal Topography 
Analysis, Invest Ophthalmol Vis. Sci., 1994; 35: 

10 2749-2757" (Reference 2), which is incorporated 
herein by reference. 

KCI is an index evaluating a similarity to 
keratoconus. KCI is found by the method 
illustrated in Fig. 4. The method of Fig. 4 is 

15 described as an expert system classifier in 
Reference 2 above. The KPI (Keratoconus 
Prediction Index) in the figure is described in 
Reference 2, which is incorporated herein by 
reference. 

20 These 18 indexes are taken as input data, 

and the probability (%) of the classification 
of each of various corneal topographies is output 
using a neural network (step 9) . Fig. 5 is a block 
diagram of the network for judging the 

25 classification of each corneal topography. In 
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this example, the corneal topographies are 
classified into known classifications (i.e., NRM 
(normal cornea), MRS (myopic refractive surgery), 
HRS (hyperopic refractive surgery), AST (corneal 
5 astigmatism), PKP (penetrating keratoplasty), 
KCS (keratoconus suspect), PMD (pellucid 
marginal degeneration), and KC (keratoconus)) 
and OTH (other known or unknown classifications) , 
i.e., 9 classifications in total. There is also 

10 provided a KSI (keratoconus severity index) 
output for evaluating the severity of KC if the 
result of the decision is KC . 

The probability of each classification of 
each of the corneal topographies is determined 

15 using a neural network (step 9) in combination 
with an expert system. The unaccounted for 
probability in the neural network process is 
represented by the OTH classification. 

The decisions of the classifications of the 

20 corneal topographies and their probabilities can 
be found as described next by entering the 18 
indexes into the neural network. 

A total of 1825 OPD-Scan® topography maps 
were acquired that encompassed the 8 

25 classifications: normal, astigmatic, 
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keratoconus suspect, keratoconus, pellucid 
marginal degeneration, penetrating keratoplasty, 
and hyperopic and myopic refractive surgery. 
These topography maps were clinically classified 
5 into one of the 8 corneal topography 
classifications. These maps were then reviewed 
and judged again by human experts in manual 
analysis of corneal topography for correctness 
of the classification and for the severity of 

10 keratoconus, if present. After review by the 
expert human corneal topography graders, 224 maps 
were chosen for neural network training and 
another 2 2 4 maps were selected at random for 
testing. The remaining maps were not considered 

15 prime examples of the 8 different classifications 
and were not used for training or testing. Each 
classification contained approximately 28 maps 
in both the training and test sets. Selection 
of training maps was based on image quality and 

20 the existence of features deemed typical for each 
classification. When KC was judged as being 
present, the human experts also assigned a KSI 
value for the purpose of training and testing 
the ability of the network to grade the severity 

25 of keratoconus. 
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Each map was F o u r i e r - f i 1 1 e r e d to a common 
resolution and 18 corneal indexes extracted for 
use as neural network inputs (SimKl, SimK2, M i n K , 
CYL, SRI, SRC, OSI, DSI, CSI, SAI, CEI, IAI, ACP, 
5 AA, SDP, CVP, KPI, & KCI) . These corneal indexes 
were previously shown to be sensitive to various 
topographical features. While 18 indexes are 
used as inputs into the total network, it was 
determined during network training that some of 

10 the individual networks that comprise the total 
network must use fewer inputs to correctly make 
a decision. Specifically, as shown in Fig. 5, 
the PKP network does not use KCI as an input and 
the KCS network does not use SIMK1, SIMK2 , MINK, 

15 SRI, SRC, ACP, I A A , and CEI as inputs (Fig. 5) . 

With the exception of the Keratoconus 
Severity Index (KSI) network which had only one 
output graded from a severity ranking of 0 (no 
sign of KC ) to 1 (maximum KC severity), each of 

20 the individual networks has two numerical outputs 
ranging between 0 (no response) to 1 (the maximum 
level of the condition seen in the training set 
examples) for each category. The first output 
is the desired classification and the second 

25 output is called "Other" , Information in the 
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first output is used for classification, while 
the information in the "Other" output is 
discarded, as indicated by the symbol "X" in Fig, 

5) . These "Other" outputs are used to make the 
5 networks train faster and more reliably. When 

the classification output value approaches 0, 
the network is indicating a total lack of 
similarity of the input variables to the input 
data originally used to train the network for 

10 that classification. When the classification 
output value approaches 1, the network is 
indicating a high degree of similarity of the 
input variables to the original training data 
for that classification. Finally, the output 

15 value is converted into a percentage, so that 
the graphical display output (as shown in Fig. 

6) ranged from 0% to 100%. This percentage 
indicates the similarity, and therefore the 
likelihood or probability of the analyzed 

20 topography pattern of the tested cornea matching 
one or more of the classifications or the 
keratoconus severity index (KSI) learned by the 
neural networks. 

Each of the individual networks was designed 

25 using a b a c k - p r o p a g a t i o n , feed-forward neural 
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network training process with a single hidden 
layer in which the network output was compared 
to the desired outcome specified in the training 
set. The single hidden layer contained the same 
5 number of neurons as the number of inputs (e.g., 
18 or fewer, depending on the category) . If the 
output error was greater than the desired 
tolerance during training, the internal 
weighting parameters of the network matrix of 

10 neurons was systematically adjusted until an 
acceptable answer was found for all training 
examples. In addition, the network was tested 
simultaneously with an independent data set after 
each training run through all training example 

15 maps. As the individual networks improved in 
the ability to correctly classify both the 
training and test set examples, the acceptable 
error tolerance in both the training and testing 
regimens was systematically set lower in order 

20 to achieve an even higher accuracy. When the 
individual network reached the maximum number 
of correct responses with the independent test 
set and no further training lead to an improvement 
with the test set, training was considered 

25 complete. This training process was repeated for 
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all 8 classification networks and for the 
keratoconus severity index network. The specific 
levels of accuracy that was achieved with the 
independent test set data are shown in the 
5 following Table 1: 



Table 1 

Category Accuracy 

KS I 9 4% 

10 KC 9 7% 

KCS 94% 

PMD 9 6% 

MRS 98% 

HRS 9 9% 

15 AST 9 6% 

N RM 9 6 % 

PKP 94 % 



The final neural network matrices that are 
20 used to define the knowledge acquired by the 
trained neural network are arrays of weighting 
coefficients. These coefficients are those 

applied at each neuron interconnection within 
the network. The mean output of each individual 
25 network can be described by a simple linear 
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equation which indicates the relative strength 
and sign of each input term toward determining 
the appropriate classification. In other, words, 
each equation is an expression of the mean sum 
5 of the weighted coefficients of each index within 
the training set. These equations are shown in 
Fig. 8 . 

Once the neural network outputs are found 
for all categories, the results are passed 

10 through an expert system threshold filter to 
obtain the final classification score. The 
threshold at which a result is displayed on the 
bar graph shown in Fig. 6 is adjustable. 
Currently the threshold is set to 10% or a 0.10 

15 output response value from each of the 
classification networks. Maps classified as 
significant for the KC category will also trigger 
the display of the KC severity in the form of 
the KSI network output, (as indicated by the 

20 control shown in Fig. 5) . If the KC output is 
not above 10%, the corresponding KSI severity 
index is not displayed. 

Fig. 6 gives an example of the result of 
the analysis of classifications of corneal 

25 topographies displayed on the display 24. The 
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indexes found from corneal curvature data and 
their values are displayed on the right half of 
the display screen. The probabilities of the 
items of the classifications of the corneal 
5 topographies are graphically represented on the 
left half of the display screen, together with 
their numerical values. Since the probabilities 
of the various topographies are graphically 
represented, it is easy to understand the 

10 contents visually. The diagnosis is also 
facilitated. KSI is given in terms of a numerical 
value when KC is more than a predetermined 
percentage, such as 50%. Since MRS, HRS, AST, 
PKP, and so on are classified together with KC 

15 and KSI regarding keratoconus, it is possible 
to judge the optimal visual aid (spectacles or 
contact lenses) for the eye. Also, the decision 
as to whether refractive surgery should be 
applied can be made precisely. 

20 Fig. 7 shows an example of a display of a 

color map of the distribution of corneal 
curvatures. When a corneal topography is 
diagnosed, the diagnosis is made more precise 
by observing the distribution of corneal 

25 curvatures at the same time. For the display of 
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the color map, the polar coordinate data 
calculated in step 7 in the flowchart of Fig . 
2 are used. The structural dependence intrinsic 
to the system has been removed. The display 
5 program can be made common to topography devices 
with different Placido ring numbers. 

In the structure according to the embodiment 
described so far, an analyzing function is 
imparted to the measuring system having the 

10 corneal topography measuring optical system. The 
cornea analysis function subsequent to step 3 
in Fig. 2 may be assigned to the computer 40. 
In this case, the corneal curvature data found 
by the a r i t hm e t i c - a n d - c o n t r o 1 unit 20 are sent 

15 to the computer 40 from the communication port 
33. The computer 40 has a function of entering 
the corneal curvature data. 

The above description of the preferred 
embodiments has been given by way of example. 

20 From the disclosure given, those skilled in the 
art will not only understand the present 
invention and its attendant advantages, but will 
also find apparent various changes and 
modifications to the structures disclosed. It 

25 is sought, therefore, to cover all such changes 
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and modifications as fall within the spirit and 
scope of the invention, as defined by the appended 
claims, and equivalents thereof. 
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